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The transition from weakly damped coherent motion to localization in the context of the spin- 
boson model has been the subject of numerous studies with distinct behavior depending on the form 
of the phonon-bath spectral density, J (M) oc w®. Sub-Ohmic (s < 1) and Ohmic (s = 1) spectral 
densities show a clear localization transition at zero temperature and zero bias, while for super- 
Ohmic (s > 1) spectral densities this transition disappears. In this work, we consider the influence 
of the phonon-bath spectral density on the noneguilibrium dynamics of a quantum dot with electron- 
phonon interactions described by the extended Holstein model. Using the reduced density matrix 
formalism combined with the multi-layer multiconfiguration time-dependent Hartree approach, we 
investigate the dynamic response, the time scales for relaxation, as well as the existence of multiple 
long-lived solutions as the system-bath coupling changes from the sub- to the super-Ohmic cases. 

Bistability is shown to diminish for increasing powers of s similar to the spin-boson case. However, 
the physical mechanism and the dependence on the model parameters such as the typical bath 
frequency ujc and the polaron shift A are rather distinct. 


Introduction 

Quantum dissipation is an omnipresent phenomenon 
in diverse physical systems, racing from quantum in- 
formatiorP and qua ntu m optics, 121 to charge transfer and 
impurity relaxation ,121^ superconducting junctions,!^ and 
more, spanning diverse energy, length and time scales. 
Describing the effects of the environment on the dynamic 
response of a sub-system requires both the development 
of theoretical and computational tools as well as the 
development of simplified models necessary to account 
for the rich system-bath dynamics and thermodynamic 
phase behavior. The minimal model required to capture 
the essential physics of quantum dissipation involves a 
two-level system coupled to a bosonic bath. Perhaps 
the most studied version is the well-known spin-boson 
model,where it is assumed that the two-level system is 
linearly coupled to a harmonic bath. The effects of the 
environment are characterized by the properties of the 
bath spectral density, assumed to have a power-law de¬ 
pendence at low frequencies, w® (s G K^), with a cutoff 
at higher frquencies determined by a characteristic fre¬ 
quency, u>c- The value of “s” classifies the nature of the 
dissipative environment, often referred to as sub-Ohmic 
for 0 < s < 1, Ohmic for s = 1 and super-Ohmic for 
s > 1. 

The dynamics, equilibrium and phase behavior of the 
two-level system is governed, amongst other factors, by 
the value of “s”. A notable transition from coherent to 
incoherent dynamics in the spin-boson model is observed 
as the Kondo parameter, rj (dimensionless strength of the 
system-bath coupling) is increased or when the bath spec¬ 
tral density changes from super-Ohmic to sub-Ohmic .SEE! 


This is followed by a localization transition at high values 
of t] in the limit Wc —>■ oo and low temperature, T —>■ 0, 
for the sub-O hmic and Ohmic cases. Such a transition 
disappears foiSEMU] § > i. This rich behavior has been 
investigated by a variety of theoretical approaches in¬ 
cluding analyticaP^HUl and approximate numerical tech¬ 
niques (for an overview, see Refs. I4l6p as well as numer¬ 
ically exact methods su ch as numerical renormalization 
group (NRG) techniques j ^ ^^ the multi-configuration‘ 
time-dependent Hartree (MCTDH) approach and its 
multilayer (ML) extension, ML-MCTDHpSHUl and path- 

integral methods 

The dynamics and steady-state properties of dissipa¬ 
tive quantum systems driven away from equilibrium have 
been the center of more recent studies, e.g. for driven 
spin-boson-type systems,^ or in the context of inelas¬ 
tic tunneling in quantum point contacts and molecular 
junctions.ESHlll xhe canonical model in the latter field 
is given by the extended Holstein model, ^ in which a 
bridge level (occupied or empty, hence two levels) is cou¬ 
pled to a bosonic bath, describing the phonons, and in 
addition to two fermionic reservoirs representing the left 
and right leads. The latter are held at a different chem¬ 
ical potentials and thus provide the source to drive the 
system away from equilibrium. Most studies focused on 
steady-state properties utilizing a variety of techniques 
to describe, e.g., Franck-Condon blockadenega¬ 
tive differential resistance,or the existence of multi¬ 
ple long-lived solutions, i.e. bistability.l^^HIZI Jn addition, 
novel numerically exa ct tec hniques uncovered interest¬ 
ing transient behavior.^^^^SSHfiH A promising approach in 
this regard is the c ombination of the reduced density ma¬ 
trix approadJsmil and the multilayer multiconfiguration 
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time-dependent Hartree method in second quantized rep¬ 
resentation (ML-MCTDH-SQR)P> used to explore the 
timescales and the dynamic “phase diagram” associated 
with the bistability.l^^lllllllllxhe most striking result re¬ 
ported by Wilner et al^foi the extended Holstein model 
with an Ohmic spectral density is that the bistability per¬ 
sists on timescales exceeding the phonon-assisted tunnel¬ 
ing time along the adiabatic potential. 

At first glance, the localization transition in the spin- 
boson model and the bistability in the nonequilibrium 
extended Holstein model may result from similar physics. 
However, these phenomena are quite different in nature 
and origin. First, localization in the spin-boson model is 
a quantum phase transition strictly at T = 0, while bista¬ 
bility is a transient phenomenon persisting over a range 
of source-drain bias voltages. Second, the former van¬ 
ishes in the adiabatic limit of a slow bath (wc —>■ 0) while 
the latter thrives in this limit. Finally, in the Ohmic case, 
localization due to a degenerate ground state and the cor¬ 
responding dependence of the steady state on the initial 
state occurs only for the symmetric spin-boson model, 
while bistability in the extended Holstein model spans a 
wide range of asymmetries 

In this work, we explore the dependence of bistability 
on the nature of the phonon spectral density using the 
combined ML-MCTDH-SQR and reduced density matrix 
approach. We cover both the the sub-Ohmic to super- 
Ohmic limits. While localization and bistability show 
opposite behavior with respect to uJc, we find similarties 
with respect to the boson spectral density. Specifically, 
as s is increased above 1 the bistability diminishes. This 
transition, however, is not as sharp as the localization 
transition. 


Model Hamiltonian and Spectral Densities 

To describe the effect of different forms of dissipation 
on nonequilibrium transport in a quantum system with 
electron-phonon interaction, we consider a generic model 
Hamiltonian describing, e.g., a quantum dot or a molec¬ 
ular junction: 

H = Hs + Hb + Vsb- (1) 

Here, Hs = edd^d is the system Hamiltonian represent¬ 
ing the electronic degrees of freedom of the quantum dot 
with creation/annihilation fermionic operators d^/d and 
energy Sd- For simplicity, we assume that the quantum 
dot is represented by a single level. The bath Hamilto¬ 
nian, Hb = Hi + Hpii, is given as a sum of electron (lead) 
and phonon baths where 

He= Ekolak ( 2 ) 

keL,R 

represents the noninteracting leads Hamiltonian with 
fermionic creation/annihilation operators and 



Figure 1: The phonon spectral density, J {iS) — , 

for various values of s for LtJc = 100cm“^and A = 1000cm~^. 


”*” 2 ^ 

j ^ ' 

represents the phonon bath with creation/annihilation 
bosonic operators b^j/bj for phonon mode a with energy 
hujj. The coupling between the system and the baths is 
given by 

VsB = X! {^kdal + tlakd^^ + dU Mj (b] + bj^j 

k^L,R j 

(4) 

where is the hopping term between the system and 
the leads and Mais the strength of the electron-phonon 
couplings to mode a. The former is determined from the 
relation 

rL,i?(e)=27r \tk\^5{e - Sk), (5) 

keL,R 

with F l,r (e) = f? ^46^ — {e — ^.l,rY used to mimic a 
tight-binding chain and ^l,r is the chemical potential 
of the left/right lead, respectively. We adopt the same 
parameters for F l,r (e) used in our recent studies 
namely, a = 0.2eV and b = leV. For this choice, F = 
0.16eV is the maximum value of rfl(£:) -|- ri(e). The 
electron-phonon couplings, M„, are determined from the 
relation: 

J (w) = TT MjS {Huj — Hujj) (6) 

3 

where we follow the notation of Caldeira and LeggettP^ 
for the phonon spectral density: 

Ji^) = (7) 

In the above equation, the dimensionless Kondo param¬ 
eter, rj = determines the overall strength of 
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the electron-phonon couplings, Wc is the characteristic 
phonon bath frequency and T (s) is the Euler Gamma 
function. For future reference, we introduce an addi- 

tional parameter A = ^ ^ f 17'^(^)> which is 

known as the reorganization energy (or polaron shift). 
Fig. shows the results for J (w) for different values of s. 

The model introduced above and variants thereof have 
been widely used to study nonequilibrium charge trans¬ 
port in nanostructures, such as, for example, semicon¬ 
ductor quan tum d ots,^ carbon nanotubes^ or molecu¬ 
lar junctions! ^° l ^^tf ^ In the latter case, the phonons may 
include, in addition to the phonons of the contacts, the 
vibrational degrees of freedom of the molecule. In all 
previous studies, however, the spectral density was lim¬ 
ited to the Ohmic case, for which s = 1. The sub-Ohmic 
(s < 1) and super-Ohmic (s > 1) limits which play an 
important role in the related spin-boson model,® have not 
been studied in transport junctions. 


Reduced Density Matrix 

The dynamics generated by the above Hamiltonian are 
rich and rather complicated to solve. Nume rically exa ct 
techniques include real time path integration^^^JsMIElHZll 
and ML-MCTDH approach.^SHHIgoth are limited to rel¬ 
atively short times and cannot describe the dynamics on 
all relevant timescales. Recently, we have proposed to 
combine a numerically exact i mpurity sol ver with a re¬ 
duced density matrix formalism .15ZE2HSMS1 Application to 
the above model for the Ohmic spectral density uncov¬ 
ered a fascinating behavior with rich dynamics on mul¬ 
tiple timescales and bistability persisting on ti mescale 
longer than the phonon-assisted tunneling times. 

Here, we adopt this approach to study the influence of 
different forms of the phonon spectral density on the dy¬ 
namic response and on the bistability. For completeness, 
we briefly review the formalism. 

The basic quantity of interest is the reduced density 
matrix, a (t), which is derived from the full density ma¬ 
trix, p(t), by the application of the projection operator 
P = pbTtb- Here, the index B refers to the bath de¬ 
grees of freedom or the “irrelevant” part of the full Hamil¬ 
tonian. cr (t) = TrBp{t) obeys a generalized quantum 
master equation, given byP^ 


initial conditions for (t( 0), an occupied and unoccupied 
dot. We assume a non-correlated initial state forpB(O) 


PB (0) = Pph (0) 0 p[(0) (g) pf (0), 
where (/3 = is the inverse temperature) 


= exp 


k^L/R 


is the initial density matrix for the leads, and 


( 10 ) 


( 11 ) 


f'ph (0) = ®xp 


-/3' 


hujj (^b]bj + i 





( 12 ) 


represents the initial density matrix of the phonon bath. 
We also consider two different initial conditions for the 


phonons, one where Sj = 0 in Eq. (12 \ corresponding to 
phonons initially equilibrated with an unoccupied dot, 
and another where Sj = 1 corresponding to phonons equi¬ 
librated to an occupied dot. More details can be found 
in Ref. 

The calculations of the memory kernel in Eq. ([^, K(t), 
is the tricky part. Formally it is given by 


R{t)=TrB{Cve-i^^*QCpB} 


(13) 


where Q = 1 — P, P = pB{0)TrB{- ■ ■ } and £ = [H,- ■ ■]. 
A more suitable form for the memory is given in terms 
of a Volterra equation of the second type, removing the 
complexity of the projected dynamics of Eq. (13l: 


K (t) = (t) — ^ (t) Cs + T / dr^ {t — t) K (r) (14) 

B Jo 

with (since the operator Cy appearing in the equation for 
$(t) and the full Hamiltonian conserve the total particle 
number, only the diagonal matrix elements <i>(t) need to 
be computed): 


ih-^a {t) = Cscr (t) +'d{t)-^ [ Stk (t) a {t - t) . (8) 
dt h Jo 

where £s = [Hs, • • • ] is the system’s Liouvillian, 

r9(t) = TrB{£ye-tQ^‘Qp(0)} (9) 

depends on the choice of initial conditions and vanishes 
for an uncorrelated initial state (which is the case consid¬ 
ered below), i.e., when p(0) = (t( 0) (g) pb{0), where cr(0) 
and p_b(0) are the system and bath initial density matri¬ 
ces, respectively, and £„ = \Vsb, ■ ■ We consider two 


$™,mm(t) = |pb {m\ ^ tkd{t)a\{t) \m) | .. 

(15) 

Here, \m) denotes the electronic state of the quantum 
dot, where m can take the values 1 or 0, corresponding 
to an occupied or an unoccupied dot, respectively. Pre¬ 
viously, we have shown that ^nn,mm(t) can be expressed 
in terms of the sum of the left and right 

currents: 


e$nn,mm(t) = -1^(0 + I^it) 


(16) 
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s=l s=l/2 



Figure 2: Transient dynamics of the average quantum dot population for the Ohmic (left set of panels) and sub-Ohmic case 
(right set of panels). The results are shown for all 4 initial conditions: Black and red curves correspond to unoccupied / 
occupied dot at phonon initial condition 5j = 0, whereas blue and green curves correspond to unoccupied / occupied a.tSj = 1, 
respectively. In all results the cutoff time used to generated the memory kernel is tc = lOOfs. 


where 

= tk{m\d{t)al{t)\m), (17) 

keL,R 

is the left/right current for an initial occupied (m = 1) 
or empty (m = 0) dot, and e is the electron charge. 

While the calculation of the memory kernel requires 
a solution of the time-dependent left and right currents, 
it typically decays on tim escales much faster than the 
density matrix itself,*^^*^ and thus is amenable to a nu¬ 
merically exact impurity solvers. For this purpose we 
adopt the ML-MCTDH approaclp 2 H 22 l calculate the 
memory up to a cutoff time tc, where is large enough 
such that the results for the reduced density matrix do 
not change for large cutoff times. When suitable, i.e., 
for small electron-phonon couplings,^ we often use the 
nonequilibrium Green’s function formalism within the 
self-consistent Born approximation to generate the mem¬ 
ory kernel for t < tc- 

Sub-Ohmic and Ohmic cases 

In Fig. 1^ we show the populations dynamics for the 
Ohmic (s = 1) and sub-Ohmic ( 5 = 5 ) cases for vari¬ 
ous values of the reorganization parameter A and the 
phonon frequency ojc- Throughout this work we fix the 
electronic parameters: Dot energy Ed = 0.25eV, ap¬ 
plied bias-voltage /i^ = —^r = O.OSeV and zero tem¬ 
perature {k-^T = 0). The results shown for the Ohmic 
case (left se t of p anels in Fig. [2 1 summarize our previ¬ 
ous findings .E3111 short, the decay of the population 
is characterized by three distinct timescales. At short 
to intermediate times, the population dynamics are gov¬ 
erned by the system-lead hybridization fi,/r(as clearly 
seen for 6j = 0, A = 1000cm“^) or by uj~^ (as clearly 


seen for Sj = 1, A = 1000cm“^). In addition to the short 
and intermediate timescales associated with the separate 
electronic and phononic degrees of freedom, the electron- 
phonon coupling introduces longer timescales related to 
the tunneling between the two charge states, as clearly 
evident for increasing values of A. On timescale longer 
than this tunneling time, the population is characterized 
by two distinct, long-lived, solutions for large values of A. 
This “bistability” vanishes as the system become less adi¬ 
abatic (increasing ojc) or for small polaron shifts (small 
A). Readers interested in a comprehensive discussion of 
the Ohmic case are encouraged to consult Ref. Ell 

The situation is similar for the sub-Ohmic case (right 
set of panels in Fig. |^. All three timescales are clearly 
observed even for s = i. For A = 1000cm“^ < Ed, the 
more stable solution is associated with the non-shifted 
configuration (Sj = 0) and the system relaxes to a single 
long-lived state regardless of the choice of initial condi¬ 
tion (Sj = 0 or 1). The dynamics at short and intermedi¬ 
ate times are governed by h/F for the non-shifted phonon 
bath (Sj = 0), since the dot energy Ed = 0.25eV is outside 
the conduction window (A/i = O.leV) and thus the sys¬ 
tem remains in the uncharged state. This is not the case 
for the shifted bath initial preparation (Sj = 1). Since 
the effective dot energy (id = Sd — 2A, Ref. EU is within 
the conduction window, we observe a rapid uncharging 
decay associated with H/T followed by a slower relaxation 
to the stable uncharged configuration along the general¬ 
ized bath mode, characterized by 1/uJc timescale. 

As A increases above Ed, the long-time behavior of the 
population depends on the choice of the initial phonon 
preparation, leading to a bistability (two distinct long 
lived solutions). The bistability is rather sensitive to 
the characteristic phonon frequency, and can lead to 
two steady states solutions in the adiabatic limit, as 
LOc 0.1^ ' As Wc increases the bistability gradually dis- 
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Figure 3: Transient dynamics of the average quantum dot population for the super-Ohmic case. The results are shown for 
all 4 initial conditions: Black and red curves correspond to unoccupied / occupied dot at phonon initial condition Sj — 0, 
whereas blue and green curves correspond to unoccupied / occupied atSj = 1, respectively. In all results the cutoff time used to 
generated the memory kernel is tc = lOOfs. Converging the results ujc = 500cm“^, A = 2000cm“^for s = 2 and ujc = 100cm“^, 
A = 2000cm“^for s = 3 is difficult, and thus we only show the direct ML-MCTDH-SQR to t = tc- 


appears. Comparing the sub-Ohmic case to the Ohmic 
case, the bistability persists over a wider range of frequen¬ 
cies and polaron shifts for s < 1. This can be explained 
qualitatively by looking at the spectral density function 
shown in Fig. In the sub-Ohmic case, the bath spec¬ 
trum is shifted to the lower frequency end as compared 
with the Ohmic case that has the same characteristic 
frequency tOc- Thus, a comparison of the dynamics gen¬ 
erated for the Ohmic case at a particular values of Wc 
should be done with the dynamics generated for the sub- 
Ohmic case (s = 1/2) for a higher uJc- 

To summarize the sub-Ohmic case, we find that the be¬ 
havior is similar to the Ohmic case. The major difference 
is observed for the long time relaxation. The sub-Ohmic 
case shows a pronounced bistability at larger values of 
LOc and a wider range of polaron shifts compared to the 
Ohmic case, resulting from the increase in the density of 
low frequency modes as s is decreased below 1. 


Super-Ohmic case 

The trends observed going from the sub-Ohmic case to 
the Ohmic case continue smoothly as s is increased above 
1. In Fig. 1^ we plots the transient population dynamics 
for the super-Ohmic case (s = 2 and s = 3) for the same 
set of frequencies and polaron shifts shown in Fig. In 
two cases, we are unable to converge the results based 
on the memory formalism with the given input and thus, 
show the direct calculations based on the ML-MCTDH- 
SQR approach generated for t < tc = lOOfs. Similar to 
sub-Ohmic case, for A = 1000cm“^ the population dy¬ 
namics at short and intermediate times is governed by 
h/T for the non-shifted phonon bath {Sj = 0) and by 
the typical phonon frequency for the shifted bath initial 


preparation {Sj = 1). The major effect associated with 
changing s is the change in the characteristic phonon fre¬ 
quency, given by sujc- As s increases above 1, the char¬ 
acteristic frequency increases leading to a rapid decay of 
the population at intermediate times, evident in Figs. 
and Concerning the long time behavior, we find that 
the tunneling dynamics are washed out as s increases 
above 1 and the bistability is limited to a narrower range 
of polaron shifts, which eventually disappears for s > 3 
for the range of frequencies studied. 

Comparing the behavior of the bistability to the local¬ 
ization transition in the spin-boson model, it is clear that 
the origin of the two phenomena is quite different. While 
in the spin-boson model localization is a quantum phase 
transition, observed strictly at zero temperature and zero 
bias, bistability is a dynamical phenomenon observed in 
the time-domain. Furthermore, while localization van¬ 
ishes for s > 1, the bistability transitions smoothly across 
the Ohmic case. Moreover, the two phenomena differ 
and show opposing behavior with respect to the depen¬ 
dence on uJc and to a large extent with respect to the 
polaron shift A. For the latter, the bistability is limited 
to a certain window of A, a window which smoothly nar¬ 
rows down as s increases. Localization in the spin-boson 
model persists above a certain value of A. 

A notable feature of the super-Ohmic case is the pro¬ 
nounced oscillations in the population compared to the 
over-damped Ohmic and sub-Ohmic behavior, as shown 
in the left panels of Fig. This is similar to the coher¬ 
ent dynamics observed for the super-Ohmic case in the 
spin-boson model. The period of oscillations increases 
as the characteristic frequency of the boson bath de¬ 
creases, which makes the convergence of the memory for¬ 
malism difficult for uJc values between 100 and 500 cm“^ 
due to the relatively long transient dynamics. As Wc in¬ 
creases, the coherence is eventually quenched because of 
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Figure 4: Left: Dot population in the super-Ohmic case for £d = 0.25eV, A = 2000cm“^ and Uc = 500cni~^. The results are 
shown for all 4 initial conditions: Black and red curves correspond to unoccupied / occupied dot at phonon initial condition 
Sj — 0, whereas blue and green curves correspond to unoccupied / occupied atSj = 1, respectively. Right: O' {a (t)} for different 
value of s for ujc = 100cm“^ (upper panel) and ujc = 500cm“^(lower panel). 


the more efficient energy exchange between the electron 
and phonon degrees of freedom, as shown in Fig.|^ 


S=2 



Figure 5: Coherent to incoherent transition in the super- 
Ohmic case s = 2. Shown are results of the dot population 
for different characteristic frequencies uJc of the phonon bath 
as indicated in the plot and for four different initial condition 
as indicated by the different colors. The color code Is the 
same as In FiglD 


solution to a (t) of the form: 

^ r (1 -b s) 

X exp [—i (1-b s) arctan (wct)]. (19) 

In the right panels in Fig.j^we plot the imaginary part of 
a (t) for two values of the frequency and for different val¬ 
ues of s. a {t) shows a transition from a smooth function 
of time at small values of s < 1 to an oscillatory function 
at larger values of s > 1. The origin of this oscillatory 
behavior is in the term exp [i (1 -|- s) arctan (wct)] and is 
correlated with the under-damped dynamics of the pop¬ 
ulation shown in the left panels of Fig. 

A more quantitative picture of the oscillatory behav¬ 
ior emerges by considering the dynamics of the reac¬ 
tion mode {Q (t)) = .^Despite the fact 

that the dynamics of the phonons have been traced out 
by considering the reduced density matrix of the system 
alone, {Q (t)) for the non-shifted case {Sj = 0) can be 
inferred using its equation of motion, resulting iiP^ 

{Q{t)) = Qi{t) + Q2{t) (20) 

with 


To understand this behavior, we propose to examine 
the phonon bath autocorrelation function, given by (at 
T = 0): 


1 f°° 

a (t) = — / duj J (oj) exp {—iujt). (18) 

“Tt Jo 


Using the spectral density of Eq. (|^, we obtain an exact 


-^5 /" Q;(T)dT (21) 

a(0) Jo 

-^5 /" 6n{T)a{t-T)dT, (22) 

a(0) Jo 

where rioo is the steady state dot population and 5n {t) = 
(Jt {t)d{t)) — nao- A similar expression not given here 
can be derived for the shifted case {Sj = 1) (see ap¬ 
pendix). The first term, Qi, on the right hand side of 
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Eq. (201 describes the relaxation of the reaction mode 
due to the coupling to the phonon bath in the absence 
of coupling to the leads. It can be calculated explicitly 
using Eq. (191, resulting in 


Qiit) 


2A ( cos(s arctan (wct)) 


shujr 


(1 +t2a;2)5 


- 1 ,(23) 


where the relations 2a(0) = 77w^r(l + s), a{T)dT = 
—iXjh and rj = 2X/{tiuJcT{s)) were used. The second 
term, Q 2 {t), includes the influence of the coupling to the 
leads on the dynamics of the reaction mode. The long¬ 
time limit of the reaction mode is determined by the first 
term in Eq. (201, i.e. 


lim {Q{t)) = lim Qi{t) = -riao 

t—>oo t—¥oo 



(24) 


In Fig. we plot the dynamics of the reaction mode 
for the same parameters used to generate the data in left 
set of panels Fig. for both the sub-Ohmic and super- 
Ohmic cases. Note that in both figures, we are unable 
to converge the results for s = 2 beyond the cutoff time 
tc = lOOfs (see the discussion above). For s < 1, the 
reaction coordinate {Q (t)) decays monotonically to its 
steady state value if it is initially in equilibrium with the 
electronic state of the dot (full black and green lines). 
The timescale governing this decay can either be purely 
phononic, arising from the term Qi(t) with an algebraic 
relaxation characteristics given by ~ (1 + 
associated with the relaxation of the electronic popula¬ 
tion 5n(t) arising from the term Q 2 {t). For s = ^ the 
decay of {Q it)) is significantly slower than that of Sn (t) 
and thus is determined by the slower phononic dynamics, 
while for s = 1 there is no clear time scale separation. In 
case of a preparation, where the phonon degrees of free¬ 
dom are initially not in equilibrium with the electronic 
state of the dot (full red and blue lines), a more complex 
transient dynamics is seen, which involves electronic and 
phononic contributions and time scales. 

In both cases (s = ^ and 1), the dot population as¬ 
sumes for longer times two distinct values depending on 
the initial conditions for the phonon bath, i.e. exhibits 
bistability, as clearly evident in Fig. The bistability 
in the dot populations leads to two solutions for {Q (t)). 
In the long-time limit, the difference between the solu¬ 
tions for {Q (t)), obtained for shifted and unshifted initial 
preparation, is related to t he co rresponding difference of 
the dot population An^ bjiSHH 


This relation between And and AQ is only valid in the 
steady state and can therefore be used as a consistency 
check if the steady state has been reached. For the results 
in Fig. [^it is fulfilled. 



Figure 6: The reaction coordinate {Q (t)) for several values of 
s and for Ed = 0.25eV, A = 2000cm“^ and uOc = 500cm“^ and 
temperature T = 0. Black - unoccupied (0) = 0 and 5j = 0; 
Red - occupied Ud (0) = 1 and Sj = 0; Blue - unoccupied 
^d (0) = 0 and = 1; Green - occupied Ud (0) = 1 and 
Sj = 1. In addition, the dashed lines depict Qi{t), the part 
of {Q (t)) which describes the relaxation of the reaction mode 
due to the coupling to the phonon bath in the absence of 
coupling to the leads (see Eq. 113). 


As s is increased above 1 the picture changes qualita¬ 
tively. First, the bistability gradually disappears as dis¬ 
cussed in details above. This implies that {Q {t —>■ oo)) 
assumes a unique value regardless of the initial prepara¬ 
tion of the phonons. Second, as clearly evident in Fig. 

{Q (t)) shows a pronounced oscillatory behavior corre¬ 
lated with the oscillations observed in the dot population 
shown in Fig. 1^ Most interestingly, for the unshifted ini¬ 
tial phonon preparation {Sj = 0), the oscillation in {Q {t)) 
lead to transient values for the reaction coordinate that 
are associated with the shifted position, and vice versa. 
This “exchange” of positions is also reflected, to a smaller 
extent, in the dot population shown in the correspond¬ 
ing panels of Fig. For larger values of s (s > 1), the 
bare phonon relaxation time scale (1 -I- is sig¬ 

nificantly faster than the electronic dynamics Sn {t), and 
thus the latter dominates the dynamics of the reaction 
mode {Q{t)) for longer times. 


Conclusions 

In this work we have studied the role of the spectral 
density of the phonon bath on the relaxation dynam¬ 
ics and bistability in the nonequilibrium extended Hol¬ 
stein model. To this end, sub-Ohmic, Ohmic and super- 
Ohmic spectral densities were considered. The results 
show a physically rich behavior, including a transition 
from incoherent to oscillatory dynamics and a disappear¬ 
ance of bistability signatures upon increase of the power 
s of the spectral density. Some features observed for the 
sub- and super-Ohmic cases can be rationalized qualita¬ 
tively by rescaling the characteristic phonon frequency. 







































Wc —>■ suJc- For example, the bistability which is exem¬ 
plified in adiabatic limit ojc —>■ 0, persists over a larger 
range of frequencies for sub-Ohmic spectral density while 
the opposite is true for the super-Ohmic case. However, 
the appearance of a slow algebraic decay for {Q it)) in 
the sub-Ohmic case as well as the oscillatory behavior of 
{Q (t)) in the super-Ohmic case cannot be explained by 
simple scaling arguments, but require a more elaborate 
analysis in terms of the bath correlation function, a (t) 
and are distinct from the Ohmic limit. 
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The equation above can be re-written in terms of the 
bath autocorrelation function given in Eq. (18 \ using the 
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1 /-oo 

= -J 7r^M2<5(a;-a;j)e-“‘da; 

1 

= - / J(w)e-*‘^‘dw 

TT Jo 

= a (t) 

and 

/■* 

y = -i / dra (r) + A, (A6) 

j J 0 

to reduce the formal solution of {Q (t)) into 


Appendix A: Reaction Coordinate Formal Solution 


Q{t) 



Sn (t) a {t 


t) dr 


In this appendix we describe the formal solution of the 
reaction coordinate Q (t) for both initial conditions (for 
simplicity we set h = 1). For the phonon Hamiltonian 



(A7) 


Hpn = (b]b, + + dUy^M, {b] + b,) , (Al) 

3 ^ ^ i 

the equation of motion for bj (t) reads 

bj (t) = —iujjbj — iMjdJd (A2) 

with the formal solution 

bj {t) = bj (0) - iMj [ Ud (r) (A3) 

Jo 

We will define the dimensionless average position of the 
phonon mode j with {xj (t)) = 


cor- 


responding reaction coordinate {Q (t)) = —^==f such 
that: 

= . -* ■ (A4) 




Using Eq. (A3) we find 


The equation above corresponds to the non-shifted 
bath initial condition i.e. {bj (0)) = 0, while for the 

shifted-bath (^{bj (0)) = , the solution reads: 


(Q(t)} = -1 


:( 0 ) 


A + 




-k 


a(0) 

I 2 

a(0) 


3 f Sn (t) a {t — t) dr 

Jo 


/o a (r)dr. 


(AS) 


In the above Sn (t) = {ud (t)) — n^o with rioo = 
limt_,.oo {ud (t)) and 3?{...} (3 {...}) are the real and imag¬ 
inary parts respectively. 

One can look at the steady state limit using the fact that 
limt_,.oo fo ('^) « (t - r) dr 0 since Sn(T) 

and a (r) approach zero for r —>■ cx) and 


poo 

/ a (r) dr = —iX 

Jo 


(A9) 


























9 


to obtain 

AQ = —AAn 

V SAUJc 

(AlO) 

where AQ and Arid are the steady state difference be¬ 
tween the shifted and the non-shifted phonon bath for 
the reaction coordinate and the dot population respec¬ 
tively. 
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